CONTINUOUS BIORTHOGONALITY OF THE 
ELLIPTIC HYPERGEOMETRIC FUNCTION 



V. P. SPIRIDONOV 



Abstract. We construct a family of continuous biorthogonal functions re- 
lated to an elliptic analogue of the Gauss liypergeometric function. The key 
tools used for that are the elliptic beta integral and the integral Bailey chain 
introduced earlier by the author. Relations to the Sklyanin algebra and elliptic 
analogues of the Faddeev modular double are discussed in detail. 
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1. Introduction 

Classical and quantum completely integrable systems serve as a rich source of 
special functions. Complexity of these systems correlates with the complexity of 
solutions of the Yang-Baxter equation - rational, trigonometric, or elliptic |22| . In 
the theory of special functions this structural hierarchy is reflected in the existence 
of the plain and g-hypergeometric functions [IJ and their elliptic generalizations 
[l9] . One of the aims of the present paper is to clarify certain questions of such a 
correspondence at the elliptic level. 

In [TTl [19] , the author introduced an elliptic analogue of the Gauss hypergeo- 
metric function V{ti, . . . ,1^) generalizing many special functions of hypergeometric 
type with "classical" properties [1]. In particular, it was shown that this func- 
tion exhibits symmetry transformations tied to the exceptional root system Ej and 
satisfies the elliptic hypergeometric equation. We start this work with showing in 
Sect. 2 that the ^-function satisfies a simple biorthogonality relation correspond- 
ing to the continuous values of a spectral parameter. We use for that the elliptic 
beta integral [16] and an integral analogue of the Bailey chains introduced in [18] . 
Our construction can be considered as an integral generalization of Rosengren's 
approach to the elliptic 6j-symbols [TT] . 
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In Sect. 3, we investigate the generalized eigenvalue problem for a finite difference 
operator of the second order D{a, b, c, q). Connection of such problems with the 
biorthogonality relations was considered in [23j . A problem analogous to ours was 
investigated earlier in ^TOl E] for the discrete spectrum. Its general solution defines 
a basis in some space of meromorphic functions. For a simple scalar product, 
we derive the dual basis defined by solutions of a similar spectral problem for 
conjugated operators. As a result, the overlap of the basis vectors with their duals 
coincides with the V^-function. In Sect. 4, we derive a new contiguous relation for 
the elliptic hypergeometric function and the biorthogonality condition associated 
with this relation. 

It turns out that the V^-function is directly connected with the Sklyanin algebra 
- one of the central objects in the quantum inverse scattering method [14l[15]. This 
relation is established in Sect. 5 on the basis of the observation due to Rains [91112] 
that the _D-operator can be represented as a linear combination of all four generators 
of the Sklyanin algebra. In this way we obtain a particular infinite dimensional 
module of this algebra with the continuous values of the Casimir operators. 

It is well known that the Sklyanin algebra can be degenerated to the quantum 
algebra Uq{sl2), which was intensively investigated during the last two decades. In 
[3], Faddeev has introduced the modular double - the quantum algebra Uq{sl2) 'S) 
Ug-i (5^2), where g is a modular transform of q (see also [B]). Originally, it was aimed 
at dealing with the well defined logarithms of quantum algebra generators, which 
can be traced to the demand of analytical uniqueness of the algebra representation 
modules. Applying the latter requirement to the results of Sect. 5, we come in 
Sect. 6 to two different elliptic analogues of the modular double defined by direct 
products of two copies of the Sklyanin algebra with the generators related to each 
other by permutation of certain parameters. An interesting fact is that some of the 
generators of these two algebras anticommute with each other. 

Finally, in Sect. 7 we discuss the Bethe ansatz for the _D-operator eigenfunctions 
and relations to the Calogero-Sutherland type models - common topics in the theory 
of integrable systems. We demonstrate that the standard eigenvalue problem for 
the D-operator can be interpreted as a difference analogue of the Heun equation. 
The latter equation is known to emerge in the one particle sector of the Inozemtsev 
integrable model [5j . Our finite difference analogue of it turns out to be related to 
the model introduced by van Diejen 2 for a restricted set of parameters. 

Before passing to the derivation of these results, we describe our notation. For 
an abelian variable z e C and a base p = e^'^*'^, Im(r) > (i.e., \p\ < 1), we define 
the infinite product 



3=0 



Elliptic theta functions with characteristics are defined by the series 



^■KiT{k+a/2f ^2Tii(k+a/2){u+b/2) 
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where the variables a and b take values or 1 and u G C. The standard Jacobi 
theta functions are defined as 

ei{u\T) = e^iu) = -Biiiu), 

02{u\t) = 02 (m) = 010 H = 0i(u + 1/2), 

03iu\T) = 03{u) = 9oo{u) = e"^/4+™0i(«, + 1/2 + t/2), 

6»4(m|t) = 9iiu) = Ooiiu) = -ze'^'^/'*+'^''"6»i(7i + t/2). 

ft is convenient to set Oaiui, . . . , Uk) ■= Oa{ui) ■ • • 9a{uk) and 9a{x ± y) := 9a{x + 
y,x-y). 

The 01 (u)-function is odd, 9i{—u) = —9i{u), and obeys the quasiperiodicity 
properties: 

01 + f ) = -01 (w), 01 (« + r) = -e-'^*^-2^™0i(u). 
The short theta function 

e{z;p) := {z;p)oo{pz~'^;p)oo, zeC*, 
is related to 0i (u) by the Jacobi triple product identity 

0i(«) = ipV8e--»(p;p)^0(e2-»;p). 
We shall need the duplication formula for 0i-function 

^ ' " + 2 ' " + 2 ' " - ^ 

and the addition formula 

6\{u±x,v ±y) — 9\{u±y,v ±x) = 9\{x ±y,u± v) (1.1) 

or ^ 

9{ux'^ jVy'^-.p) — 9{uy^ ,vx^;p) = —9{xy'^,uv^;p), (1.2) 

where 9{ux^;p) := 9{ux,ux^^;p) := 9{ux; p)9{ux^^ : p) . 
The standard elliptic gamma function has the form 



= n - i^'i' i«i < 1- (1-3) 

We set 

rp,q(il5 • • • ) ^fc) -fp,g(^l) * * * ^p.,qi^k)i 

Tp^q{U^) ■■= ^pAizKAtz'^)^ rp,g(2±2) ._ Tp^^(^z'^)Tp,g{z-^). 
This function satisfies two finite difference equations of the first order 

^P,qiiz) = 6{z;p)rp^q{z), Tp^qipz) = 9{z; q)Tp^g{z) , 
the reflection equation 

^p,g{z)^p,giP<l/z) = 1' 

the duplication formula 

^pA^') = rp,«(-^, -z, q'/\ -qy\py\ -py\ {pqy/'z, -{pqf'^z), 

and the limiting relation limp^o ^p,q{z) = 9)00- 

With the help of three pairwise incommensurate parameters a;i,2,3 G C, we define 
three base variables 

27ri^ 2ni^ 

q = e , p = e '"^ , r = e "i, 
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and their particular modular transforms 

q — e "^1, p — e ^3 ^ r = e ^3 _ 

The modified elliptic gamma function, which remains well defined for UJ1/UJ2 > 
(i.e., \q\ = 1), has the form 

G(w;c.) =r,,,(e2"^)r,v(re-'"-) = e-"^^(")r(e-'"^;f,p), (1.4) 

where \p\, \r\ < 1 and 

It is fixed uniquely by the following three equations 

fiu + LO,) = 0(e2--/-. . /(^ + = 0(e2-»/-i ; r)f{u), 

/(u + c.3) = e-"*^^'^("^"^'"^)/(^i), (1.5) 

B2,2(U;iJl,UJ2) = h h h - 

UJ1UJ2 LOl LO2 bLU2 bUJl 2 

and the normalization /(X]L=i ^"i/^) ~ ^- reflection equation has the form 
G{a;uj)G{b;cu) ^ 1, a + b = Y:L,ujk. A more detailed description of the elliptic 
gamma functions can be found in reviews |13j or [19j and in the references therein. 

2. Continuous biorthogonal functions 
We investigate elliptic hypergeometric integrals of the form 



j('")(^i,...,w6)^^ / , (2.1) 

2m+6 



where T is the unit circle and \tj\ < 1. For m = 0, we obtain the elliptic beta 
integral [Hj 

/(°)(ii,...,t6) = n ^p^^i^jtk). 

l<j<fe<6 

Using this formula, it is easy to verify the following recursion relation by mere 
change of the order of integrations 



Xk/ T^pr. 1^ '[ti, . . . ,t2m+i,£mW,ernW )— {,2.2) 

Jt rp,g(w±^) W 

where Em = \/Ylt™^2m+5 ^k/pq- It leads to the m-tuple integral representation for 
/("i) similar to the Euler type integral representation for the general plain hyperge- 
ometric function m+iFm, i-e. Z'™^ can be interpreted as an elliptic analogue of the 
m+ifm-function. Recurrence (|2.2p is a special realization of an integral analogue 
of the Bailey chains discovered in [TB] (the Bailey chains technique is well known 
as a simplest tool for proving the Roger s-Ramanuj an type identities [T|). 



n2m+5<fe<;<2m+8 ^p,q{tktl) 
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For m — 0, the integrand on the right-hand side of relation (12. 2|) contains the 
elhptic beta integral, computation of which yields [17] 

Vih,...,ts)^ n ^p..q{tkti,tk+4ti+i)Vieti,...,eU,e~H5,...,e-Hs), (2.3) 

l<fc<i<4 

where V{t) = /'-^-'(i.) is an elliptic analogue of the Gauss hypergeometric function 
and e — ^ tc^tQtjt^ / pq — \/pq/tiht^. This is the key reflection extending the 
obvious S's-group of symmetries in parameters of V{t) to the Weyl group for the 
exceptional root system E-^. 

For m = 1, we obtain the equality 

t(2)u i \ Tl7<k<l<10^P-<l(^l'^l-) 

^(,^l,...,^loj = , 2x 

r^lU^) V ft,,...,te,e,^,.r^-^)-, (2.4) 



where ei = y^trtstgtio/pq. After fixing parameters t^ti — tets = pq, the integral 
in the left hand side gets reduced to the elliptic beta integral and, so, it can be 
computed explicitly. This yields the relation 

Y[ ^p.A'tktl) Yl ^p,q{tkt9,tktlo) = n7<fc<i^l0 ) 

i<fc<(<4 fc=i Tp^,{tsho}Tp,g(e^) 

XK — - — j_ V{ti, . . . ,U,pq/t7,pq/t8,eiw,eiw ) — , 

where tii2i3^4^9iio = PQ- We denote now tg = s^,tio — s^~^ and obtain 

l<fe<i<4 j=l ^P.^l^li 

Jt rp,,(w±2) 

X y(ti, . . . ,t4,pq/tr,pq/ts,eiw,eiw~^) — , (2.5) 



where tit2ty,tiS^ — PI and £i = s\/tjt^Jpq. After the replacements 

ti — c, 12 — a, i3 = s \ — -X, 14 — 3 \ — - X , tj ~ a, tg = o, 

V cd y cd 



the condition tit2t^t4s'^ = pq is satisfied automatically and ei — s^J ab/pq. Intro- 
ducing the basis vectors 

cj^iw; a, 6|C; s) = Tp^.i-sa^^ , sb^^ , \/!|^^^*)' ^^'^^ 
we can rewrite the key relation (|2.5p as 

0(a;; c, s) = K / R{c,d,a,b;x,w\s)4>{w,a,b\S^;s) — , (2-7) 
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where, after the substitution s — y/pq/p, 



R{c, d, a, b; x, w\s) 



d ' 



X V I c, d, aI —X, ij —X 
cd V cd 



-1 pq pq ab 



w, 



-w 



p 



(2.8) 



Using the transformation (|2.3p . we can bring this i?- function to the form 
i?(c, d, a, 6; X, w\s) 



X y sc, srf, 



-X, e\ 



[pq -I pq pq o.b 

' as ' fes ' y pq 



The kernel R does not depend on the variable ^. Relation (|2.7p looks therefore 
like a rotation in some functional space spanned by the basis vectors 0, which is 
performed by the "matrix" R with the continuous indices x and w. 

We can denote ^ — e'^, w — e^^ , and a; = e"'' and let x, and take real values 
only. Then, because of the symmetries 

(l){w;a,b\£,;s) = (/)(w~^ ; a, s) = (l){w;a,b\C^;s), 

we see that this function depends actually on the variables y := cos6 and cosx- 
Similarly, the function 

R{a, b, c, d; x, w\s) = R{a, b, c, d; x~^,w\s) = R{a, b, c, d; x, w~^\s) 

depends on y and cos tp. Using the elementary relation 



/(cos^)^ 



f{cose)de 



1 vi-y 



we can write 



x;c,d\S^; s) — 2iK / R{c,d,a,b;x,e'^^\s)(j){e^^ ;a,b\£_; s) 



dy 



y 



(2.9) 



Taking the limits c a and d ~f b and denoting v — cos ip, we obtain the relation 



27rv/r 



hm i?(c, d, a, 5; e^^ e^'^) = ^ , " 6{v ~ y) 

c^a,d^b [p;p)oa(q]q)oa 

in the distributional sense. Formally, this equality can be rewritten as 



V a,b,. 



pq pq I ab 
a ^ b ^ 



(p;p)oo{q;q) 



S{v-y). (2.10) 



Using twice equality (j2.7p with different choices of parameters on the right-hand 
side, we obtain immediately the self-reproducing property for the kernel 

f div 
K I R{a,b,c,d;x,w\s)R{c,d,e, f;w,z\s) — = R{a,b,e, f;x, z\s) (2-11) 

. T W 
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and the biorthogonality relation 
fi 



R{a, b, c, d; e^^,^'\s)R{c, d, a, b; e'^'l^) = ^^^-^ <5(z; ~ 

where v — cos ip and u' — cos can be considered as continuous spectral variables. 
Substituting the explicit expressions for the i?-functions, we rewrite the biorthog- 
onality relation as 



1 



V\a,b, jP-e^^, J^e-'"^, ^, ^, W-e^^ W-e^^ 



X V I c,d,jAe''',J^,e- 
V ca 



a6 

pq I ah 



pq pq I cd 
c d \l p 



I cd 
P 



-tip 



dy 



a 



ab p cd p ±2i(p\ Vl - v' 



^'"^ ^ p ^ ab^ p ' cd' 



(2ik) 



■ (5(w — v'). 



(2.12) 



3. A SCALAR PRODUCT AND BIORTHOGONALITY OF THE BASIS VECTORS 

It is necessary to establish now the Hilbert space content of the construction 
described above. It is desirable to connect the F-function with a scalar product of 
some basis vectors of this space. 

Multiplying equality (|2.7p by 0(z;pq/e,pq//|^;s^^) and integrating it over ^ 
with the measure k Jj^pq{£,^^)d£_/£_, we obtain on the left-hand side 



Jt 



(z;f ,f |$;s-i)0(x;c,d|e;3)de 



rp,,(C±^) 



{sti,...,st4,S ^t5,...,S ^s) , (3.1) 



where ti = c, <2 ~ f^, ^3,4 = a; s ^y/pq/cd and — pq/e,tQ — pq/fih.} 



z^^ S\fcflpq. On the right-hand side, we find 

K / R(c,d,a,b;x,w\s)V (st\, . . . ,st'A,s~^t5, . . . ,s^^tg) — , 

where t[ = a,t2 — b,t'^ ^ = w^^s^^ y/pq/ab. If we apply iJy-transformation ()2.3|) to 
T^-functions on both sides of the equality, we obtain relation (|2.1ip , which shows 
the self-consistency of the consideration. 

Taking the limits e ^ c, f ^ d, and applying equality (|2.10p to the right-hand 
side of (|3.1|1 . we find the biorthogonality relation 



j<p' . m m\p. „-i 



f |e;s-i)0(e'^;c,d|e;s)dC 



rp,,(e±2) 

P<1 cd ±2^y 



a/i — (5(w — v'). 



(3.2) 



The latter result inspires introduction of the following inner product. For two 
functions of z € C with the property '4'(z) — xjj^z"^), we define 



{X{z),%l){z)) = K 



X{z)ip{z) dz 



(3.3) 
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From the reflection equation for the eUiptic gamma function, we find for 2; e T and 
real (p, (7)-variables 

1 2 

i.e. our measure is positive definite. Similarly, for p = q* we obtain a positive 
measure as well 

l^p,q[Z ) 

We introduce now the finite difference operator ^ [TT] 

9{az,bz,cz,dz;p) 0{az-'^ ,bz-'^ , cz'^ , dz-^;p) -^^/^ 

where T^q^'^f{z) — f{q^^^^z) is a g-shift operator. The operator conjugated to Dr 
is defined by the equality 

{x,DBia,b,c,d;p;q)ijj) = {D*jj{a,b,c,d;p;q)x,tp)- 

Assuming that functions "tpiz) and x{^) do not have singularities in the annulus 
< \z\ < |(7|-i/2, we find 

^ . , ^ n M^/^ pq^^^ pq^^^ pq^'"^ 

Dj^(a,b,c,d;p]q) = Dr , — — , , —j—]p]q 

\ a b c d 

Clearly, D*^ = Dr. For further considerations, it is convenient to scale one of the 
parameters by p and define D(a, b, c, q) := —aDp>{pa, b, c, d;p; q), so that 

6l(az,6z,cz,dz;p) g(az~\bz-\cz~\d2:~^p) 

(3.4) 

and 

D (a, 6, c, d-p- q) = -^^D _ ,j . 

We look now for solutions of the generalized eigenvalue problem for ZJ-operators 
with the fixed product p :— abed, 

D{a, b, c, d;p; q)f{z) = XD{a, b, c', d';p; q)f{z), (3.5) 

where cd = c'd' . We parametrize the spectral variable A by an elliptic function of 
the second order 

(wy^ c/d, W\fdjc\ p) 



A = Xiw) = , , , 

e{w./djdf,w./dljd-,p) 

where it; S C is an arbitrary parameter, and indicate in f{z) dependence on all 
parameters 

/(z) -.^ f{z-w-q^'\,q''\p), 

where the parameters a and 6 are scaled by q^/^ for a later convenience. Equation 
p.Sp can be rewritten then as 

f{q^/^z-w]q^/^a,q^''^b-p) _ _^^^d{z^ , az-'\bz^^;p) 



f{q~^^^z; w; q^/'^a, q^/'^b] p) 0{z^^, az, bz\p) 

9{cz^^, dz^^ ,W\J d /d' , Wyjd' /c';p) — 6{c'z~^, d'z~^ , w^/c/d, Wyjd/c;p) 
9{cz, dz, W\fd~[d! , w^J d' /c';p) — 9{c'z, d'z, W\fcjd, w^J d/c]p) 
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which, after appUcation of the addition formula for (?(z;p)-functions, reduces to 



f{q-^/^z; w; q^/^a, q^/%; p) 0{az, bz, V^w^z;p) 

A particular solution of this equation valid for jgl < 1 has the form 



(3.6) 



f{z]W]a,b]p) ^Tp^q\—z ,—z '-l/— ^ ^ I- (3-7) 

The general solution at [gj < 1 is obtained from (|3.7p after multiplying it by arbi- 
trary function Lp{z) satisfying the equality (p{qz) — ip{z). The case |g| — 1 will be 
considered below, but the regime |g| > 1 will be omitted for brevity, since it turns 
out to be equivalent to the case |g| < 1 after a change of parameters. 

Our particular solution depends only on three parameters a, 6, and p, i.e. two 
variables, say c and c', drop out being plain gauge parameters. A special choice of 
parameters reduces f{z; w; a, 6; p) to products of theta functions 

f Uq\ Ar;a,b;q-^\ =l[e{q^az^;p) [] e{q^bz^;p), k,NeN, (3.8) 

which coincide with the intertwining vectors of Takebe [H] and the basis vectors 
used by Rosengren in [TTl [T2] . 

Since the function f{z]W]a,b]p) is symmetric in p and q, we can demand it satisfy 
one more generalized eigenvalue problem obtained from (jS.lOp by permutation of p 
and q: 

D{a, b, c, d; q;p)f{z; w;p^/^a,p^/^b; p) 

0{Wy^,Wy^;q) 1/2 1/2, VQOA 

= —, rrfif ^— — D{a,b,c ,d;q;p)f[z;w;p ' a,p ' b;p)(3.9) 

ff{wy'cJ/a', w^d' jc; q) 

where cd = c'd' . For incommensurate p and q (which we assume throughout the 
paper), this requirement strongly restricts the space of solutions of our eigenvalue 
problem since the equations (p{qz) — ^p{pz) — ip{z) are satisfied only by a constant 
function. Our basis vectors /(z; w; a, b; p) become thus fixed uniquely up to a 
constant multiplier which we set equal to 1 for concreteness. 

Function p.7p satisfies another type of the generalized eigenvalue problem 

D{a,b, c, d;p; q)f{z;w, q^^^a, q^^^b; p) — d^^d{w^ ^/d/c, cd; p) f {z; w; a, b; p) . 

(3.10) 

Indeed, dividing this relation by /(z; w; a, b), we obtain on the left-hand side 

0{az, bz, cz, dz;p) /{q^^'^z; w; q^^^a, q^/'^b; p) 
z9{z'^;p) f{z;w;a,b;p) 

6{az^^,bz^^,cz^^,dz^^;p) f{q^^^^z; w; q^^'^a, q^^^b; p) 



(z 



-2. 



p) f{z;w;a,b;p) 



After substitution of the explicit form of /(z; w] a, b; p), this expression reduces to 



9{cz , dz , z ab/ p; p) — 9{c ^z,d ^z, w'^z^J p/ab]p) 
'79{z^) ■ 

Addition formula p.2p with u = z\fcd, x — \fcjd, v — z^/abjp, y = w yields an 
expression required for the right-hand side in p.lOp . Equality (|3.5p is easily deduced 
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from (I3.10|) by noticing that f{z;w; a, b; p) does not change after replacements of c 
and d by c' and d' with the constraint cd = c'd' . 
The following commutation relation [TU] 

D{a, b, c\ d';p; q)D{q^/^a, q^'\ q-^'\, q-^'^d-p; q) 

= D{a, b, c, d-p; q)D{q'/^a, q^'\ q-^'^c' , q^^'^d'-p- q), (3.11) 

where cd — c'd' , is a simple consequence of p.lOp . Indeed, the action of both sides 
of this equality on the function /(z; w] qa, qb; p) yields the same result. The straight- 
forward verification of the identity is also possible (coincidence of the terms in front 
of the operators TJ*^^ is evident, and the scalar terms coincide after application of 
the addition formula for theta functions). 

Consider the generalized eigenvalue problem for the conjugated operators 

D*{a,b,c,d;p;q)g{z) ^ \{v)D* {a,b,c' , d';p; q)g{z), 

where cd — c'd' and X{v) is a new spectral variable. Evidently, this equation is 
obtained from (13. Sp after replacements a ^ a = pq^/^/a, b ^ b = pq^^^ /b, c-^ c = 
q^^^/c, d ^ d — q^^'^/d, and w — > u. We find thus its particular solution 

g{z) := g{z; v; a, b; p) = Tp^q ^az=^, 6z^, \J^'"^^^^ ■ (3-12) 

The biorthogonality relation of interest emerges as follows 

= {{D*(a,b,c,d;p;q) - \{v)D* {a,b,c' ,d' ;p; q)) g{z;v;a,b; p), 
f{z;w;q^^^a,q^^^b;p)) = {g{z;v; a,b; p), 

{D{a, b, c, d;p: q) - X{v)D{a, b, c', d';p; q)) f{z; w; q^^'^a, q^''^;p)) 

= (^1 - (3(^5 ^! P)^D{a, b, c, d;p; q)f{z; w; q^^'^a, q^/'^b; p)) 

= {^^ TTt) d^^^i"^^ '/Wc,cd]p){g{z]V]a,b] p), f{z]W]a,b] p)\?,.li) 



\{w) 

For \{v) ^ A(w), we find thus {g{z\ v\ a, 6; p), /(z; w] a, b; p)) — 0, which essentially 
coincides with relation (|2.10p . Restrictions on the parameters needed for conjuga- 
tion of the operator D without crossing singularities are \p\ < |a|, \b\ < \q\^^^ and 
\p/abq\ < \v\^ < \abq/pl \ab/pq\ < \w\^ < \pq/ab\. 

To conclude, solutions of the generalized eigenvalue problems for the ZJ-operator 
and its conjugate define a basis rotated by the i?-function and its dual. An overlap 
of these bases defines the ^-function: 

g{z;v;a,b;e)f{z;w;c,d;e) dz ( ± pq pq /— ±1 .01 ^^ 

\a,b,J—v^,—,—,\ —w^ \ , (3.14) 




Tp^z^'^) z \ ' ' V a& 'c'd 

where V{. . . , ax^, ■ ■ ■) = V(. . . , ax, ax~^, . . .). 

We can rewrite also the key relation (|2.7p in a compact form 



Tp^az ^ ,l3x ^ ) = n I r{a, P,"f,6; z,x;t,w)rp^q{'yt ^ ,6w ^ )- 

Jt ""^ 

r(a,/?,7,(S;z,a;;i,w) = - — — ^ — (3x^,az^, —t"^, — 



where a/? = jS. 
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4. Elliptic hypergeometric equation 

We recall the derivation of the elliptic hypergeometric equation together with 
its general solution [121 120] ■ Repeating transformation (|2.3p with the parameters 
e^a, £^4, £"^^5, £"^^6, playing the role of ti,. .. ,1^, and permuting t^,t4^ with 1^,1^ 
in the result, we find 

4 

V{t:} = n ^pAtjtk+i) V{VT/h, Vt/U, Vu/h, Vu/ts), (4.1) 

j,k=l 

where T — tit2^3i4 and U = t^tQtjts- Equating right-hand sides of (|2.3p and (|4.ip . 
we obtain the third symmetry transformation 

Vis)^ II rp,,(s,Sfc)1/(VM/si,...,VM/s8). (4.2) 

i<i<'i-'<8 

All these relations represent different Weyl group transformations for the excep- 
tional root system Ej in the parameter space. 

If we set V ^ tQ,x = tj^y = z, and u = in the addition formula p.2p . then it 
becomes equivalent to the equation 

= t69{tstf;p)A{z,ti, . . . ,te,qtr,ts), 

where A{z,t) — Yit^i ^p.q{tkZ^) /^p,q{z^'^) is V^-function's kernel. Integration over 
z along an appropriate contour C yields 

M {ts4;p) V{qte) - te0 (tstfip) V{qh) = t^O {tet^;p) V{qts), (4.3) 

where V{qtj) denote V^- functions with the parameters tj replaced with qtj and 
the balancing condition 11^=1 ~ P^l- Substituting transformation (|4.2p in this 
equation, we obtain 

hO Y[Hkh;p) V{qtr, qts) - ^ ^ (*7ifc;p) V{qte, qts) 

= te9 l[0{t8tk;p)V{qte,qt7), (4.4) 

where 11^=1 *i = P^- 

We permute now parameters and tg in the latter equation, change the param- 
eters ^5,^7 — > t^/q,tj/q, and exclude the function V{qtQ,q~^tY) with the help of 
equality (|4.3p where the parameter is replaced with t^/q. The resulting relation 
looks like 



0{t7tt/q;p) , 
0{Utt-p) 



V{q-'h,qts)~Vit) 



0{tj/qte,t5ts/q,qts/t7;p) -A- 0{t5tj/q;p) 



0{tstf,qte/t5;p) fj^ 0{t7t,/q;p) 

4 



9{qts/t7,tj/t5;p) Yj 0{tQtj;p) 



0{h/t&,qhlh]p) ^Jj 0it7tj/q;p) 
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Multiplying it by the function 

.f. . . e{ts/qt(i,tets,h/ts;p) yr e{trtj/q;p) 

and symmetrizing in parameters and ig yield the elliptic hypergeometric equation 

A{h, ...,tr, h;p) (^^^7^^^^Viq-%, qts) - Vit)) + [tj ^ tg) + VH) = 0, 

(4.7) 

which was derived in [19j in a slightly different way. 
We set tr = cx~^ , ts = cx and denote 

T>{ti, ...,ts)^ a{ti, . . . ,tj,ts){Ta:^q - 1) + a{ti, . . .,ts,tr){T~g - 1) + K{ti, . . .,ts), 

(4.8) 

where T^.qfix) = f{qx) and 



e{ts/qte,tets,te/ts:p) 



5 



&[t8 h,h qts,p) 



k=l 



K{ti ,...,ts)= e{t7ts/q; p) n ^(ieifc; p) ■ 



k=l 

Equation (|4.7p may be rewritten then as a linear second order g-difference equa- 
tion V{t)U{t) = 0, where U{t) = V{t)/Tp,q{t'jt^ ,ts,t^). Other linearly independent 
solutions of this equation can be obtained by replacing tj — > ptj, for some fixed 
j = 1, . . . , 5, or c ^ pc, OT X ^ px, etc [19]. 

It is not difficult to verify the operator identity 

_ e{t3/t'„t,/t'^;p) _ h 0itytltyt'i;p) 

where in the primed set of parameters t' variables ^3 and are replaced by arbi- 
trary and t'^ with the constraint t^ti = and similarly for the double primed 
parameters. From the elliptic hypergeometric equation 'D{t)U{t) — 0, we find that 
the ratio 

V{t')U{t) _ V{t")U{t) 

e{u/t'^M/A\p) ^ e{ult'iM/t'l\p) 

is independent on the auxiliary parameter tg (or t'^. To find its explicit form, we 
multiply contiguous relation (|4.5p by A{tl) and symmetrize it in the parameters 
and ig. After some tedious computations, we obtained 

vmm_ o{^.tMt.,^-^:p) 

r- — t^tQ y r- l^.iUj 

*2. *3. • 1 f) { *3t7 tats t4t7 t4t8 qte . \ 

V 9 ' 9 ' 9 ' 9 ' ts ' 



X 



n (^[—■^P)U{1 h,qt,) 



( gfe t3t7 fill- T) 

V ts ' 9 ' g 



' ihte,Ute;p) 



U{t). 
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The elliptic hypergeometric equation easily follows from this relation. Indeed, it is 
sufficient to replace tg by t'^ in ()4.10p . to subtract resulting equation from (|4.10|) . 
and to use operator identity (14. 9p . 
In the new parameter notation 

e/c = — , K = 1, . . . ,5, £6 = cieP , £8 = —, £7 = — 
ctk h q 

with the balancing condition Hfe^i £fe — P^Q^ > the operator D takes the form 

V^iei, . . . ,£8) A{x)iT.,^g - 1) + A(x-i)(r-i - 1) + z^, (4.11) 

where the coefficients 

obey the explicit ^e-symmetry in £1, . . . , Eg. Evidently, the equation 'Dx{e)'ip{x; e) = 
has a solution 7/1(2;; £) = U(ti, . . . , te, cx, cx~^), where 



tfe = ,A; = 1,...,5, <6 = — , c= 5 — . 

C£fc £8 

The operator Vxie) has the following conjugate with respect to the inner product 

8 

X'*(£i, . . . ,£8) = a(x)X':r(ei, . . . ,e8)a(a;)"\ a(a;) = J| rp^^(£fcX=^), (4.13) 

k=l 

provided the functions -0 and (j) do not have singularities in the annulus \q\ < \x\ < 
A solution of the equation 'D*{e)x{x; £) = is therefore given by the function 

x{x;e) = a{x)U{ti, . . . ,te,cx'^,cx). 

We denote as e' the set of parameters e where £3,4 are replaced with £3 4 such that 
£3£4 = £3£4 (or, equivalcntly, ^3^4 — > tjj 4, ^3^4 = ^3^4). The evident equality 

= {K(M:)x{x;s'),P{x)i^ix:^)) = {x{x;s')Me')P{x),pix;s)), (4.14) 

where P{x) is some function, generates an infinite family of biorthogonality re- 
lations. For instance, if /? = 1, then the function x(x;£') is orthogonal to the 
function on the right-hand side of (|4.10p for 6'(£3/£3, £4/£3) 7^ 0. This complicated 
biorthogonality is clearly different from (|2.12|) . However, it is plausible that for 
some particular /3(x) relation (|4.14p yields equality (|2.12p . 

5. Connection to the Sklyanin algebra 



The Sklyanin algebra [111 [15] is generated by four operators So, ... ,83 satisfying 
relations 

SaSp — SpSa — i{SoS^ + SySo), 

SoSa — SaSo = iJpj{SpSj + SjSp), (5.1) 

where the triple (a, /?, 7) represents any cyclic permutation of (1, 2, 3). 

Sklyanin has found the following realization of this algebra in the space of theta 
functions of the even order 2N: 
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where e^''^" are the shift operators, e^"^^^ f{u) — f{u ± 77), and 

Sq{u) = 6*1(77, 2u), si{u) = e2iv, 2u) ^ ei{7] + 1/2, 2u+ 1/2), 

S2(u) = ie^ir,, 2u) = ^e"("/2+^+2")0l(r7 + 1/2 + r/2, 2u + 1/2 + r/2), 

S3(u) = 6'4(?7, 2u) = _e"*(^/2+7,+2«)^^(^ ^ ^/2, 2u + t/2). 

The structure constants Jap in this case are parametrized as 

They can be represented in the form Jaf) = {Jp — Ja)/J-y, where 

_ ^2(277)02(0) _ 03(277)^3(0) 04(277)^4(0) 

0^(77) ' 02(77) ' ^IW • 

Verification of the defining relations is rather comphcated and requires multiple 
use of the Jacobi identity [12] 

20i(6-B/2) = 0i(fe) + 0i(6+l/2)+e"*("+^'(0i(fe + r/2)-0i(6+l/2 + r/2)), (5.2) 

where S = 61 + 62 + 63 + ^4 and 0i(6) = nt=i &iih)- 

There are two Casimir operators for this algebra. The first one is 

a=0 a=0 J // 

-0a+i(2u - 25, -27/ -2g- 277)) + (m ^ -7/) = 402(2.9 + 77). 

The terms proportional to e^^'jSu cancel due to the Jacobi identity with the choice 
61 = 62 = 77, 63 = 27i — 2g, ~ 2u — 2g + 277, so that B ~ A{u — g + rf) and 
9i{b~ B /2) — 0. The scalar terms combine to A0\{2g-^ri) due to the same identity 
with the choice 61 = 62 = 77, = 2u — 2g, 64 = -~2u — 2g — 2r] with B = —Ag. 
The second Casimir operator has the form 

a—1 a—1 ^ ^ 

-0a+i(27i - 2g, -27i - 2g - 277)) + (u ^ -u) = 40i(25, 2g + 2f]). 

The terms proportional to e^^^^" cancel due to the Jacobi identity with the choice 
= 0,62 = 277,63 = 274-25,64 = -27^ + 25-277, so that B = and0i(6-B/2) = 0. 
The scalar terms combine to 40i(25, 25+277) due to the same identity with the choice 
bi = 0, 62 = 277, bs — 2u — 2g, b^ — —2u — 2g — 2r] with B ~ —Ag- We stress that 
all the described operator relations are satisfied for arbitrary continuous parameter 
5 (i.e., the choice 5 = N'rj/2 of [15] corresponds to a special finite dimensional 
representation in the space of theta functions) . 
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As noticed by Rains [TOl [T2] , a linear combination of the Sklyanin algebra gen- 
erators Sa is equivalent to the Z)-operator. Indeed, it is related to the combination 

n'=i^iK+«4+25) n'=i^iK+«4+25+i) 

2A(ai,a2,a3,a4) — So — — Si 



71 



2^ 

n?=i0i(a,+a4_+25 + i±^) 

l+r ^ 



■ +2a4+2g-r,) llj = l"iV"-J ' i ' 2 ^ c 

ei('?+^) 



01(^+i) 

where we assume that ai + a2 + 03 + 04 = —4(7. Substituting explicit expressions 
for Sa and using the Jacobi identity with bi = ai + a4 + 2g, 62 = 02 + 04 + 2g, ^3 = 
as + a4 + 2(7, 64 = ~2{u — g) (so that B = 2a4 + 4g — 2u), we find 



rij^i^iK nj=i0i(% -") 



A(ai,a2,a3,a4) = ^ ' fl^ o ^ 

Conversely, any Sa can be viewed as a special case of the A-operator: 

S2 - Xe-01 (v+'-^]A('-±^-g,'-^-gJ-^-g ' 



4 -"4 4 4 

^3=Xe-0i(.+ 0A(l-,,-I-,,i + I-,,4-I-, 

where x = W^^^ / {p]p)oo, P = e^'^*'^. 

In the multiplicative system of notation 

(a,5,c,(i) e2"'^l'=•^'^ p := a5c(i 6"^"^ z ;= e277™^ g.^^i^ir,^ 
the operator A can be rewritten as 

A(ai, a2, 03, a^) = (ip^^^ip;p) 00^ e'*"^L'(a, 6, c, d;p; q), 

where the operator D coincides with (|3.4p . In the paper [TU] it was indicated that 
equality (j3.1ip comprises all commutation relations of the Sklyanin algebra. The 
solution of our generalized eigenvalue problem (|3.10p , fixed in (|3.7p , defines thus a 
special infinite dimensional module for the Sklyanin algebra with continuous values 
of the Casimir operators. 

Relation (|3.10p describes the action of the /^-operator on this module f{z; w; a, b; p) 
with a special choice of parameters. To know the action of the Sklyanin algebra 
generators in general, we need to consider 

D{a, b, c, d; q;p)f{z; w; q^^'^e, q^^'^h; p) (5.3) 

9{az, . . . ,dz, —w^z;p) 

'' -{z^ z^^) I f{z;w;e,h;p), 



6{z^;p) \ z^9{ez,hz;p) 



16 



V. P. SPIRIDONOV 



where p = ahcd. We apply now a known theta functions identity (see Lemma 6.4 
in[n]) 

n n+2 

j=i fe=i 



ai---a„ ^ n;=i,5:^j^(aj7ai;P) 



where nj=i % 11^=1 ^fe = 1- Namely, we choose n — i, bi — a, . . . ,b4^ — d^h^ — 
e^^jfli = h^^, 0,2.3 — u)^^ \J eh/ p and find after some simplifications on the right 
hand side of (15.31) 



ph jJ= f{z; q^^^'^w] qe, h; p) + {w ^ w^^) 

\ h^e{hJfw,w^;p) 

, e{h-^a,...,h-H,h-^e-^-p) \ 

Fixing parameters ai,2,3 in a different admissible way yields hidden permutational 
symmetries of this expression. In particular, it is symmetric in e and h. 

The I?-operator maps thus basis vectors to their linear combinations with shifted 
values of parameters e, /i, and w. Using relation (|2.7p with an appropriate choice 
of parameters, we can replace /(z; q~^/'^w] qe, h; p) by the action of some integral 
operator on the vector f(z;x;qe,qh]p) with integration over x. This provides a 
realization of the D-operator as an integral operator acting on basis vectors with 
the parameters qe and qh. 

In [12j . Rosengren proved the Sklyanin conjecture about the reproducing kernel 
associated with the discrete basis p.Sp and the Sklyanin invariant measure [H] . We 
use a different measure p.3p . i.e. relation (|3.14|) does not represent a generalization 
of the similar overlap product for the elliptic 6j-symbols in [T^] . 

6. Elliptic analogues of the Faddeev modular double 

Now we would like to discuss the uniqueness of the basis vectors /(z; w; a, 6; p). 
As we saw above, they are fixed up to a constant multiplier by two generalized 
eigenvalue problems (j3.5[) and (j3.9p . The first equation uses linear combinations of 
the standard Sklyanin algebra generators 

Sa = *'""^7#? {0a+i{2u - 2g|r)e''^" - 0a+i{~2u - 2g\T)e'^^-) . 
6i(2u\t) 

The second equation requires introduction of another copy of this algebra, obtained 
from the first one by permutation of q and p (or of r and 2r]) : 

Sa = {ea+ii2u-2g\2rj)ei'>- - 9a+i{-2u - 2g\2rj)e-i^-) . 

For that, we should have |g| < 1 (i.e., Im(?7) > 0). The structure constants J^p are 
also obtained from Jq,^ by permutation of r and 2r]. Evidently, the S'a-operators are 
not well defined in the limit Im(r) —^ +oo, i.e. such a way of fixing representation 
modules exists only at the elliptic level. 
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The two sets of operators Sa and Sa satisfy the foUowing simple cross-commutation 
relations 

SaSb^SbSa, a,be {0,3} or a, fee {1,2}, 

SaSb^-SbSa, a e {0,3}, foe {1,2} or a e {1, 2}, 6 e {0, 3}. 



Algebra (|5.ip has many automorphisms [15j . For real or purely imaginary rj 
(i.e., q* = q^^), the transformation S* = —Sa is an involution for imaginary r (i.e., 
p* — p). Because the parameter 77 should lie in the upper half plane for the Sa- 
operators, such an involution exists for both algebras only for purely imaginary rj. 
There are specific involutions existing only for our direct product of two Sklyanin 
algebras. For instance, at rj* — — t/2 (i.e., q* = p) the conditions Im(r), Im(r/) > 
are preserved and we obtain S* ^ —Sa and J*p — Ja/3- 

There is a second possibility for fixing solutions of equation (|3.5p uniquely. For 
its description, we use additive notation and substitute t = lu^/uj2, 2r] — uji/uj2, 



z — e ^2 , a,o,c,a = e , w — e , p = e 

and f{z; w; a, b; p) h{u; v; a, /3; a) in (|3.6p . which yields 



h{u ■ 



■,v;a + ■ 



,f3 + 



{az ^,bz ^,w^z ^\/cd;p) 



h{u -^■,v;a + ^,P+^;a) 0{az, bz, w±zV^;p) 

A solution of this equation valid for |p| < 1 and Ig] < 1 has the form 

G{a/2 + P/2 ~a/2±v±u; iv) 



h{u] V] a, /3; a) 



G{a ±u, (3 ±u;u!) 



(6.1) 



(6.2) 



where G{u;uj) is the modified elliptic gamma function. It is defined up to the 
multiplication by an arbitrary function 'p{u) such that ip(u + ui) = ^(u). The 
gamma function G{u]uj) is symmetric in toi and uj2- It remains a well defined 
meromorphic function even for iOi/uj2 > (i.e., \q\ = 1 or 77 is real). Demanding 
that function (|6.2p solves simultaneously equation (16. 1|) and its partner obtained 
by permuting lui and lu2 for uji/ll!2 > 0, we fix (p{u) — constant. 

In terms of the Sklyanin algebras, we have now two copies of it with the second 
algebra being obtained from the first one by permutation of oji with uj2- Such a 
transformation is equivalent to changes rj — > 1/4?/, r — > t/2?7, and u — > u/2r] in Sa, 
which yields 



9a+l (4^ 2rj) 



2ri 



u~ g 



277 



1] 



49„ 



We have now the following cross-commutation relations for Sa and Sb- 

SaSb = SbSa, a, fee {0,1} or a, fee {2,3}, 

SaSb^-SbSa, a e {0,1}, fee {2,3} or a e {2, 3}, fe e {0, 1}. 

Consider the degeneration of these Sklyanin algebra generators in the limit 
Im(w3/a;i), Im(w3/w2) ^ +00 (i.e., r,p — > 0), which is permissible now. It is 
not difficult to see that for p ^ expressions p~^^^So,i, p^^^{S2 + iSs), and 
— iSs) have well defined limits to some finite difference operators [1]. 
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Simplifying the latter operators by shifting u — > u + c and taking the limit e := 
g^27ric _^ obtain three operators E,F, and k: 

,1/4 _ g-1/4 - + ' .^n gl/4+g-l/4 -^-fc ' 

lim pi/8^±i^ ^ ^_2.»(g-„a„ _ ^. (gi/2 _ ,-i/2)j7^ 

where g = e^'^*''. The operators i?, -F, k form the quantum algebra Uq{sl2) with the 
defining relations 

P - 

kE=^q^/^Ek, kF = q-^/^Fk, EF - FE = —j-^ — . 

The ^a-operators degenerate in a similar way to 



g2,, _ e 



e " (e 1 _ g „ +2"") 
^ = il I 

forming the ?7q-i (s^2)-algebra 

fc2 - r 



fei; q-^/^Ek, kF = EF - FE ■ 



g-i/2_gi/2' 



with (J = e^"'/''. We obtain thus the Faddeev modular double [3j. 

In [3] , it was demanded that generators of the copies of quantum algebras com- 
mute with each other. This demand is satisfied by the explicit realization of the 
modular double used in [3]. In our case, not all the generators of two algebras 
commute with each other: 

kF = -Fk, kE = -Ek, kF = ~Fk, kE = -Ek, (6.3) 

which seems to contradict this demand. However, introducing the operator K = k^ 
and rewriting defining relations of t/g(s^2)-alg6bra as 

KE = qEK, KF = q-^FK, EF - FE ^ ^/2~_^-i/2 ' 



we obtain the desired commutativity of generators E, F, K with E, F, K, as in 
[B] (note that q and q of the paper [B] equal to our parameters g^/^ and q~^^^, 
respectively). At the elliptic level, such a quadratic change is not natural; one has 
to work with the operators Sa and Sa some of which anticommute with each other. 

For \q\ = 1, algebra Ugish) has the involution E* = -E,F* = -F,K* = K. 
As noticed by Faddeev [3], for q* = q, the modular double has another involution 
permuting algebras in the pair: E* = —E, F* — ~F, K* = K. For our direct 
product of the Sklyanin algebras, in the first case we have S* = —Sa and S* — —Sa, 
provided J*^ = Jap and J*^ = Jap (or, equivalently, \q\ — 1 and p* — p)- An 
elliptic analogue of the Faddeev involution has the form S** — —Sa, and it exists 
under the conditions J*^ — Jap- These constraints for the structure constants are 
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satisfied, if \2rj\ = 1 (i.e., 77* ~ l/Arj or q* ~ q) and r* ~ — t/2?7 (i.e., p* = r). 
Equivalently, 2r/ = e"^, G (-tt, tt), and r = ipe'"^/^, p > 0. 

Summarizing our consideration, we conclude, that the second way of fixing func- 
tional freedom in the definition of basis vectors provides an elliptic generalization 
of the Faddeev modular doubling principle. The first way of fixing the functional 
space does not have an evident trigonometric limit, i.e. at the elliptic level there 
are two different "modular doubles" . 



7. Bethe eigenfunctions and relation TO THE Heun equation 

As shown by Sklyanin, all operators Sa preserve the space of even theta functions 
of order 2N, if 

The operator A(ai, . . . ,04) obeys thus the same property (its multivariate gener- 
alization to the root system BCn proposed by Rains [5] obeys similar property as 
well). Therefore one can try to diagonalize it in the space of such functions: 

N 

A{ai,...,a4Mu) = Ei;{u), tp{u) ^ 9i{u±u^), (7.1) 

m— 1 

'ip{—u) = Tpiu). This equation has the following explicit form 

^ [[0,{u±u^ + r;) + [[9,{u±u^-^) 

N 4 



m=l k=l 

If we divide this relation by tp{u), then we obtain on the left-hand side two elliptic 
functions of u of the order 27V -I- 4 and a constant on the right-hand side. In order 
to verify this elliptic functions identity it suffices therefore to check its validity for 
2N + 5 different values of u in the fundamental parallelogram of periods for some 
particular choice of the spectral parameter E. 

First, it is easily verified for u = 0, 1/2, r/2, (1 + t)/2, when 9i{2u) = 0. Second, 
the validity of equation (|7.ip at u = ztuj yields N Bethe ansatz equations for the 
roots Ml, . . . , Wat: 

Oijak + Mm) _ A Oijum ±Un-r]) _ 

Actually, it suffices to satisfy any iV — 1 equations among (|7.2p , since elliptic func- 
tions cannot have only one simple pole in the fundamental region. As a (2iV-|-5)-th 
value of u for which equation (j7.ip is valid, we take u ^ ai for some I which fixes 
the parameter E as 

„ _ rifc^i + tt 6li(a; ±Un + r]) 

01 (2a,) JJ-^ Oi{ai±Un) ' 

Apparently, this expression should not depend on the choice oi I — 1, 2, 3, 4, which 
is rather easy to check for 0. 
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One can analyze in a similar way a more general Bethe ansatz 4'{u) — e~'^ Ilfc^i 
Uk) for some constant 7 determined from the eigenvalue problem A?/' — Eip. 
There are solutions different from (|7.ip . For instance, when the operator A is 
reduced to 5*0, then for N = 2n the roots Uk can be chosen in such a way that 
ip{u) — e'''" Jl^^j^ 0i(2(u — Uk)), and the Bethe equations have solutions with a 
non-trivial parameter 7 (see review ;13 and references therein). All these spe- 
cial eigenfunctions are expected to be related to some finite difference operators 
commuting with A and forming some algebraic relations with it. 

We consider now the limit 77 — > 0, corresponding to the transition from finite 
differences to differentiations. For that we perform a gauge transformation 

A(ai, ...,a4)= xp^^{u)A{ai, a4)p{u) 
^ p{u + v) 6*1 (g + u) _ pju-rj) Oija ~ u) 
^ p{u) 9i{2u) ^ ^ p(u) 0i(2u) 

where x is some constant to be specified later on. Choosing p(u) as a solution of 
the equation 

, 1 6Ii(2m) , , 
p{u -f r/) --^-f^, — 7P(w), 
X t^i(a + u) 

we find 

A(ar, . . . , a.) = e^^" - X^^^^^ft^-^e-^^. 

t^i(2u, 2u — 277) 

The potential 

B\(u — d,u ^ a — T]) 
ei{2u,2u-2r]) 

is an elliptic function of u of the eighth order with the periods 1 and t. We choose 
now 

1 T l + T 

ai = ao7], a2 = - + airj, a-^ = - + a2ri, 0,4 = h 01377, 

and consider the limit 77 — > 0. Fixing 



X ■ 



jpl/8g7ri(a2-a.'i)') 



we find after some computations A(ai, . . . , ai) — 2 — rfL + 0(77'^), where 

4=0 

and uj[) — O^uJi — 1,102 = T,ujj, = uji+uj2, and the function V{u) = —d? log 0i{u)/du^ 
coincides up to an additive constant with the Weierstrass function with the periods 
1 and r. 

The eigenvalue problem for the L-operator is known to lead to the Heun equation 
with the general set of parameters. The equation determining the spectrum of the 
A-operator represents thus a finite difference analog of the Heun equation. We call 
therefore A(ai, . . . , 04) (or D{a, b, c, d;p; q)) the difference Heun operator. 

Since the Heun equation is related to the one particle Hamiltonian of the In- 
ozemtsev model [5] (an elliptic Calogero-Sutherland type model for the i3C„-root 
system). Therefore is it natural to expect that the /^-operator has a similar re- 
lation to relativistic quantum multiparticle systems. Indeed, this is true for the 
model proposed by van Diejen [5] and investigated in detail by Komori and Hikami 
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[5. In the presence of the balancing condition 11^=1 — P^Q^^ the corresponding 
one particle Hamiltonian takes the form 

H = Ai.){T^,-l)+Ai.-'){T-l-l), M-) = ^fi^^^- (7.4) 
Imposing the additional constraints 

and replacing £5 with pes, we find Jl^^s = 1 ^^'^ 

H = -e^^ {D{e5, . . .,es;p;q'^) - ejr^6'(ese6, esey, eses;?)) . 

A more detailed consideration of relations between this model and the ZJ-operator 
lies beyond the scope of the present paper. We see thus that the difference Heun 
operator is a rather universal object: it emerges in the Calogero-Sutherland type 
models, it generates the Sklyanin algebra and a family of biorthogonal functions 
related to the elliptic analogue of the Gauss hypergeometric function. 

The standard eigenvalue problem for Hamiltonian (|7.4p with a special choice of 
parameters and for a special eigenvalue coincides with the elliptic hypergeometric 
equation [5D]. Therefore, a particular eigenfunction of the ZJ-opcrator is related to 
the F-function. To find it, we substitute ti^2 — ig^/^c"^, = ±(p(j)^/^c~^ in 
and find that under these conditions 

a{ti, . . .,tj,t8) + a(ti, . . = K{ti, . . .,ts), 

t& f c c pq 2 



V{ti, ...,t8) = D cte, — , — , —;p;q 

c \ te qte ct^ 



This means that the function 



V ( ^ u U cx-^ 



where ts^e = pc^i is a zero mode of the resulting D-operator. However, this zero 
mode is easily computable from scratch 

r 2(^x''^ ^x"^) 

ip2[x) - 



Meromorphic solutions of the equation Dip(x) = are defined up to a multipli- 
cation by an arbitrary function satisfying equality ip{q^x) = (p{x), i.e. the ratio 
ipi{x)/ip2ix) is an elliptic function of x. 

The Heun equation is the general differential equation of the second order with 
four regular singular points. Description of its general solution in the closed form 
is one of the important open mathematical problems. From our consideration we 
see that solutions of the generalized eigenvalue problem for the difference Heun 
operator defines the basis of functions, the connection problem for which is solved 
by the i?-function ([^ . 

The modular parameter r can be considered as a time variable for the Hamilton 
dynamics generated by the classical analogue of operator (|7.3p . Then the equations 
of motion are known to be equivalent to the Painleve-VI equation [S] . In a similar 
way, we can take the classical analogue of the D-operator and interpret it as a 
non-stationary Hamitonian with the time variable r. For a special dependence of 
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the parameters aj on t, the corresponding equations of motion would define some 
particular g-deformation of the Painleve-VI function. 

Biorthogonality (|2.10p is easily generalizable to a similar relation for the entire 
hierarchy of integrals (|2.ip . Therefore, it looks possible to generalize the construc- 
tion described in this paper to I*^™^ with m > 1. It is natural to conjecture also 
that there exist analogues of the basis vectors f{z] w] a, 5; p) for the root system 
BCn, which can be built with the help of the generalized eigenvalue problem for 
the Rains operator ,9j. Their determination is a challenging problem for future 
investigations. 
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